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^ Abstract 

fvj We propose a high efficiency high fidelity measurement of the ground state spin of a single NV 

^ center in diamond, using the effects of cavity quantum electrodynamics. The scheme we propose 

^__j is based in the one dimensional atom or Purcell regime, removing the need for high Q cavities 

I— I that are challenging to fabricate. The ground state of the NV center consists of three spin levels 

'j^ ^^{m=o) ^iid '^^{m=±i) (the ±1 states are near degenerate in zero field). These two states can 

3 undergo transitions to the excited (^E) state, with an energy difference of ~ 6 — 10 /^eV between 



:3 



the two. By choosing the correct Q factor, this small detuning between the two transitions results 
in a dramatic change in the intensity of reflected light. We show the change in reflected intensity 



K* can allow us to read out the ground state spin using a low intensity laser with an error rate of 

'^ ' Rs 7 X 10^^, when realistic cavity and experimental parameters are considered. Since very low levels 

^^ of light are used to probe the state of the spin we limit the number of florescence cycles, thereby 

^^ limiting the non spin preserving transitions through the intermediate singlet state "^A. 
O 
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Addressing single spins is an important route to quantum computation [T]. The long 
decoherence times of spins such as trapped atoms |21 E], ions[l], or charged quantum dots [3], 
make them ideal candidates for storing and processing quantum information. There are 
many schemes for using internal spin states in all these architectures [U Ej [7], resulting in 
the demonstration of fundamental quantum logic gates[Hl El [THj- NV centers in diamond 
have long decoherence times even at room temperature [TT] making them another promising 
candidate for performing quantum information tasks. Several experiments have shown the 
manipulation of the ground state spin of a diamond NV center using optically detected 
magnetic resonance techniques (0DMR)[TT1[T2]. This has further led to the coherent control 
of single ^^C nuclear spins and quantum logic operations [13], HI]. The main problem in using 
ODMR is that the detection step involves observing fluorescence cycles from the NV center, 
which has a probability of destroying the spin memory. Since the energy level transitions 
of the NV center are not polarization sensitive, we cannot use Faraday rotations to perform 
quantum non demolition measurements of the spin state as was shown for charged quantum 
dotsfini US]- The scheme we propose here is similar to the ODMR scheme, however by the 
introduction of a low Q cavity we vastly reduce the number of photons required to probe 
the spin state, therefore keeping the disturbance of the ground state spin to a minimum and 
not destroying the spin memory. 

If we consider the energy level structure of the NV center in figure 1, the ground state 
is a spin triplet split by 2.88 GHz due to spin-spin interactions [IH]. The excited state is a 
triplet split by spin-spin interactions, but with the further addition of spin-orbit coupling |2U]. 
Recent experimental evidence [T7] has uncovered this excited state structure (figure 1). The 
net effect of spin-spin and spin-orbit interactions is to create a detuning ^ lA GHz(6 /ieV) 

1(^=0) and the ^^(^=+1) 

l(m=-l) ^^ '~ '~^" '~^^"' 

state spin of the defect. The energy level structure is not simply a ground and excited 
triplet state, there also exists an intermediate singlet state ^A. There is a probability that 
the ^E state can decay to this state, with different rates depending on the spin state. For 
the ^Em=±i states (transitions 6,7) both theoretical predictions and experimental results 
suggest the decay rate is around 0.4 x l/r[20l [211 122], where r is the spontaneous emission 
lifetime (~ 13 ns). For the ^-Em=o state (transition 5) theoretically the rate of decay to 
the singlet should be zero[2U], however experimental observations have shown the rate to be 



between the transition from the ^^(^=0) and the ^74(^=^1) to ^E state or ^ 2.5 GHz(10 /ieV) 
for the ^Arm=-i) It is exactly this detuning that we plan to exploit to measure the ground 
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FIG. 1: Energy level diagram of the NV center in diamond showing the experimentally determined 
ground and excited state splitting|17| [T8] . The defect has zero phonon line at 637 nm, with width 
of order MHz at low temperatures [19] 



^ 10"^ X l/r[21j. Since the ^A singlet state decays preferentially to the ^Em=o state [20| [23] 
(transition 8), then it is clear from the rates above that broadband excitation leads to spin 
polarization in the spin zero ground state P^. Since transition 8 is non radiative then there 
will be a dark period in the fluorescence when it becomes populated, and as the decay rate 
from ^Em=±i to the singlet state is much larger than from ^Em=Oi the change in intensity 
measures the spin state[21|. Clearly using fluorescence intensity to detect the spin state has 
a probability to flip the spin, therefore it would seem necessary for a scheme to suppress 
this. However spin flip transitions are essential to initialize the system. Thus a compromise 
is required between the perfectly cyclic spin preserving transitions required for readout, and 
the A type spin flip transition required for initialization. 

We consider the structure in flgure 2, which can be modeled as a single sided cavity, 
where k is the cavity decay rate (side leakage), r] is the coupling of the cavity to external 
modes, g is the NV center cavity coupling rate and 7 is the NV dipole decay rate. We can 
write down the Heisenberg equations of motion for this structure as|25]: 
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FIG. 2: Schematic diagram of an NV center embedded in a photonic crystal cavity with cavity 
decay rate k coupled to a photonic crystal waveguide at a rate 77. 
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where uja and ujc are the atomic transition (cr_) and the intracavity photon (annihilation 
operator d) frequencies respectively, a^ represents a Pauli Z operator on the atomic state 
and measures the population inversion. If we now combine this with the input output 
relation for this cavity: 



hn - bout = ^f^^a, (4) 

then we can find the refiection coefficient for light input into the cavity via 6^: 
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where we have set a^ = — 1 as is appropriate for the weak excitation limit. At low temper- 
ature the zero phonon linewidth is 0.1 /ieV[T9], we set g = 0.03 meV as appropriate for a 
cavity mode volume of 0.02 /im^, where the NV has an oscillator strength of ~ 0.2 given a 
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FIG. 3: Plots showing the reflection spectra against detuning where uJa = ^c for (a) k = 75 /ieV, 
(b)K = 7.5 ;ueV,(c)K = 0.75 /ueV,(d)K = 0.075 /ueV, and i] = 50k for ah of the above 



13 ns lifetime. It is desirable for the cavity to be critically coupled to the input output so 
we will set r/ to be 50 times faster than k. 

Figure 3 shows the effect of varying k and ri on the reflection coefficient. When k and r] 
are very low we are in the strong coupling regime with g » k, r/, 7, here we can clearly see 
the two Rabi split dressed atom-cavity states in figure 3(d), and all of the light resonant with 
the cavity mode is refiected. As we increase k and 77 we are no longer able to resolve the two 
states and cross over into the one dimensional atom regime rj + k > g > '-/, where in figure 
3(b) there is a small peak in refiectivity on resonance, a result of quantum interference |26]. 
In this one dimensional atom or Purcell regime the damping of the atomic transition (zero 
phonon line) plays an increasingly dissipative role as a larger proportion of the radiative 
decay is into non cavity modes. The result of this is a dip in reflectance clearly visible in 
figure 3(a). 

In figure 4 we have plotted the refiectance at the cavity-atom resonance against rj as it is 
the dominant decay channel for the cavity mode. The amount of refiected light drops to zero 
at a value of 77 = 4(^^/7 which corresponds to the transition from the one dimensional atom 
to the weak coupling regime. Thus at this point all of the light absorbed by the NV center 
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FIG. 4: Calculated reflectance against cavity damping at uJc = ^a and o^c — u; = 0, where we have 
ignored k as r] = 50k. The relevant atom-cavity coupling regimes are labeled, the strong coupling 
regime occurs when rj is very low and |r(LL')| ~ 1. The crossover from the Purcell regime to the 
weak coupling regime occurs at a value of 36meV as predicted. 

is emitted into non cavity modes. After this turn over point we are in the weak coupling 
regime, where the NV center has progressively less effect on the dynamics of the system as 
it so weakly coupled. It is near this transition region that we wish to operate where the 
narrow feature in the reflection spectrum caused by the zero phonon linewidth dominates 
(figure 5). 

If we consider figure 5, then if we set the cavity to be resonant with the ^A(m=o) to excited 
state transition then the reflected intensity for resonant exitation becomes: 



kHlm=0 



'-/{k, — 1]) + Ag^ 



0.001 



(6) 



7(fi; + rj) + Ag'^ 

Almost none of the input light will be reflected when the NV center is in the spin m = 
ground state. However if it is in the spin m = +l{m = —1) state then the 7 /ieV(10 /xeV) 
detuning means that the NV center is effectively uncoupled giving a reflected intensity of: 
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Nearly all of the input light will be reflected. This contrast in intensity (|r(c(;)|m=±i ~ 
|r(u;)|m=o) can be easily detected. What makes this result significant is that the curve in 
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FIG. 5: Plot showing how intensity contrast (|r(ti; 



lm=±l 



r cj 



^o) varies with the ratio of r/ to 



K. Inset: The reflection spectra against detuning where Wa = ^c for k = 750 ^eV, and r\ = 50k. 



figure 5 corresponds to a total Q factor Qtot = uj/{k, -\-ri) ^ 55. Since we have set rj to be 50 
times greater than k, this means the photonic crystal cavity before coupling needs to have 
Q = uj/k, ^ 3000. This is much lower than the cavity Q factor of 300000 that would be 
required to have the cavity linewidth narrow enough to resolve the two transitions, which 
when coupled to a waveguide in the same way as here would need to exceed 10^. It is possible 
to further reduce the requirements on the cavity Q factor by reducing the the ratio of rj to 
K. However the result of this is a reduction in the intensity contrast between the two spin 
states as a larger proportion of the light confined in the cavity leaks out of the side. The 
intensity contrast which measures the spin is not infiuenced by total Q factor, the optimal 
value being Qtot ~ 55, the contrast is only infiuenced by the ratio of rj to k,. It is desirable 
to have this contrast at a maximum in order to minimize errors in state identification. 

There are several benefits to this scheme. The first is the obvious increase in collection 
efficiency of the photons, making low intensity measurements possible. Since less photons 
are required to probe the spin state there are less fiuorescence cycles therefore a reduced 
probability of a spin fiip transition. Additionally as the cavity is resonant with the m = 
transition by probing with narrow band light then we never excite the spin ±1 transitions 
which have a higher probability to spin-fiip, hence the system is optimized for spin preserving 
transitions. However if we pump with a broad band laser source we can easily spin polarize 
the ground state to initialize the system. Since we are in the low Q regime then the Purcell 



factor is small, Fp ^ 4 for a system with the parameters listed above, thus the rate of 
spontaneous emission(SE) into the cavity is not significantly modified. If we were operating 
in the high Q or strong coupling regime then the SE rate into the cavity would be much larger, 
and the decay rate to the -^A singlet state would remain unmodified. Therefore with the 
probability of a spin flip transition greatly reduced, the system would not be simultaneously 
optimized for readout and initialization. 

In order to make the scheme experimentally relevant the limitations of current detector 
technology must be included. If we consider an overall detection efficiency of 33% then if 
we input 60 photons we can expect to detect 20 with unit reflectivity. If the spin is in the 
m = ±1 state (|r(co')p ^ 92%) it is reasonable to expect 18 photons to be detected. If the 
spin is in the m = state (|r(c<j)p < 1%) we may expect 1 photon to be detected. If we 
set a detection threshold of 6 photons the error in the measurement can then be found from 
the probability of detecting > 6 photons when we expect 1 and the probability of detecting 
< 6 photons when we expect 18, giving an error rate of ~ 1.5 x 10^'^ (assuming poissonian 
distribution). Standard silicon avalanche photo diodes have a dead time of 50 ns which 
means that it will take 3 /iS to carry out a measurement with 18 detected photons (running 
at one third of the detector saturation count rate). Since the longest observed spin coherence 
time of an NV center is 600 /is[27] this introduces a further error rate of 5.5 x 10"^. There are 
also errors associated with saturation of the NV center. However as the detector dead time 
is much larger than the modified spontaneous emission lifetime then these are negligible. 
Finally there is also an error associated with decay from the ^Em=o to the ^Em=±i state via 
the ^A singlet state which for 60 photons is < 10"'^. Thus the total error rate is ~ 7 x 10~^. 

Simulations of photonic crystals in diamond have shown Q factors larger than 10^ for mode 
volumes around 0.02 /im^, or larger than 10^ for mode volumes around 0.008 /im^|28j. These 
values are significantly more demanding than those required for this scheme, particularly 
the Q factor. Experimental evidence suggests that the actual Q factors will be much lower 
than those simulated. Cavities fabricated by Wang et.al.[29j, showed more than a factor of 
10 shortfall in the experimental Q factor compared to the simulated, attributed to defects 
in the nanocrystalline structure. Nevertheless their measured Q factor of 585 would allow 
a ratio of ?7 ~ IOk, in order to have an overall Q factor of 55. This would result in a 65% 
contrast between the two spin states, increasing the error rate to ^ 2 x 10~^. Theoretical 
considerations of the absorption in nanocrystalline diamond have predicted a reduction in 
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Q factor from a value of 66300 to a value of around 1350 for a cavity of mode volume 0.02 
/xm^lSn]. For our purposes this would result in a contrast of 85%, where the error rate would 
be ~ 1 X 10^^. So the scheme is clearly robust and can cope with experimental imperfections. 
The main difficulty with this scheme, which is true for all schemes, is the positioning of the 
NV center at the field maximum. If the precision is poor then this can have a detrimental 
effect on the coupling rate g. This in turn reduces the intensity contrast between the two spin 
states, which is sensitive to the value of g compared to the zero phonon linewidth 7. There 
is promise that ion implantation in single crystal diamond could hold the key to fabricating 
suitable devices[ni], the precision of implantation is currently on the nanometer scale[2IlE2]- 
The use of single crystal diamond would dramatically reduce the absorption losses caused 
by defects, so experimental Q factors should be closer to the theoretical predictions. 

In conclusion we have proposed an efficient low error measurement of the ground state 
spin of an NV center. For the realistic parameters proposed here we can achieve error rates 
of around 7 x 10"'^. The setup can easily switch between initialization and readout by 
switching from a broad to narrow band laser source. Low error readout requires modest Q 
factors, and even with current limitations of photonic crystal cavities the error rate could 
be as low as 2 x 10~^. Work needs to be done on the design and fabrication of photonic 
crystal cavities coupled to waveguides, particularly in single crystal diamond to minimize 
absorption losses. We also note that we can measure the spin with a single photon with 
92% fidelity (assuming ideal detection), where fidelity is simply the contrast between the 
two spin states. Hence with some modifications the ideas here could be used to remotely 
entangle two spatially separated NV centers embedded in cavities, which is a subject for 
further study. 

We acknowledge support from the UK EPSRC (QIP IRC), and the European Commission 
under projects IST-015848-QAP and IST-034368 EQUIND and Nanoscience ERA project 
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